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Abstract 

An equation of state(EOS) of nuclear matter with explicit inclusion of a 
spin-isospin dependent force is constructed from a finite range, momentum 
and density dependent effective interaction. This EOS is found to be in 
good agreement with those obtained from more sophisticated models for 
unpolarised nuclear matter. Introducing spin degrees of freedom, it is found 
that at density about 2.5 times the density of normal nuclear matter the 
neutron matter undergoes a ferromagnetic transition. The maximum mass 
and the radius of the neutron star agree favourably with the observations. 
Since finding quark matter rather than spin polarised nuclear matter at the 
core of neutron stars is more probable, the proposed EOS is also applied 
to the study of hybrid stars. It is found using the bag model picture that 
one can in principle describe both the mass and size as well as the surface 
magnetic field of hybrid stars satisfactorily. 



PACS Number(s): 97.60.Jd, 21.65.+f 
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1. Introduction 

Since the pioneering works of Baade and ZwickyfllJ and Oppenheimer 
and Volkoff[f2[ about half a century ago, and the identification of pulsars as 
rotating neutron stars around the year 19680, several studies have been 
made to unravel the mystery of the structure of the neutron stars. It is 
a well accepted fact that the density inside a neutron star varies from the 
surface to the core by about 15 orders of magnitude. Understanding the 
structure of such complex objects requires an accurate knowledge of the 
equation of state(EOS) of neutron star matter in the different density re- 
gions. The extremely low density domain and the subnuclear region can be 
well described by the EOS given by Feynmann-Metropolis- Teller (FMT)[|J 
and Baym-Pethick-Sutherland (BPS)[|5j respectively. The EOS of the dense 
nuclear matter is still riddled with some uncertainties. A consistent EOS 
should describe the compressional properties @ of nuclei near the ground 
state and also of the hot and dense nuclear material that is created in en- 
ergetic nuclear collisions 0. It should also answer such important questions 
as whether a star at the later stage of its life explodes or not|| and how 
the neutron star is born. Renewed interest in the rapid cooling of neutron 
stars H [KJ by the direct URCA process has also an immediate bearing on 
the nuclear EOS. 

In the non-relativistic framework, EOS of nuclear matter have been con- 
structed with phenomenological effective interactions or from realistic inter- 
actions with different degrees of sophistication[ll-16]. Most of these EOS can 
explain the properties of neutron stars like their mass, size, moment of inertia 
etc. well within the observational limits. However, in regard to the magnetic 
properties of neutron stars, till now no acceptable explanation exists for the 
origin of the rather large magnetic field (~ 10 12 G) at the surface. A recent 
analysis [0 of binary millisecond pulsars suggests that a permanent compo- 
nent of this magnetic field could exist, sustained by a spontaneous magne- 
tised phase inside the neutron star. Attempts have been made to explain the 
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presence of the magnetic field by means of a ferromagnetic transition. 

In the framework of Hartree-Fock theory, employing hard and also soft 



core potentials, Pfarr[|18| does not get such a ferromagnetic transition. A 



similar conclusion is reached by Forseth and OstgaardJT9| who made the 
calculations in the lowest order constrained variational method of Pandhari- 
pande using soft core potentials; with a hard core potential a transition to 
the ferromagnetic state was seen to occur though at ~ 30 p Q , where p Q is the 
normal nuclear matter density. In a relativistic a + u Hartree-Fock approach, 
a ferromagnetic transition is also predicted by Niembro et al|2(| at too high a 
density. However, in an improved model [BlJ with inclusion of 7r and p mesons 
in addition to a + uj, a ferromagnetic transition is seen to occur at a compar- 
atively much lower density, p ~ 3.5p D , but the incompressibility of normal 
nuclear matter is found to be too high (~ 450 MeV). It would therefore be 
interesting to know whether a ferromagnetic phase transition is possible at 
a density realisable in neutron star matter with an EOS with firmer grounds 
in the experimental realities of normal nuclear matter and finite nuclei. This 
is the primary motivation of this work. 

There is a strong possibility that at the core densities of neutron stars, 
there is a phase transition from nuclear matter to quark matter [22-25]. About 
a decade ago, Witten conjectured f26[ that the strange quark matter(SQM) 
might be the absolute ground state of hadronic matter i.e., the mass energy 
per baryon may be less than 930 MeV. If this is true, then the possibility that 
the pulsars are rotating strange quark stars may not be ruled out. Even if 
SQM is not the absolute ground state (i.e. at densities less than the hadron- 
quark transition density the nuclear matter is energetically favoured than 
SQM), one may still find hybrid stars having quark cores with nucleon en- 
velops. Since our understanding of the confinement /deconfinement of quarks 
is far from complete, all the aforesaid possibilities are only speculative. One 
should also keep in mind that the theoretical framework used in general to 
study the phase transition is phenomenological and simplistic in nature. 

In view of all the above, we would like to investigate in this paper the 
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following. Firstly, whether there exists an EOS which consistently describes 
the nuclear matter and finite nuclear properties, as well as predict a ferro- 
magnetic transition at a density realisable in the interior of neutron stars. 
Secondly, whether the same EOS which predicts a ferromagnetic transition 
permits a hadron-quark phase transition and thereby the formation of a hy- 
brid star. And finally using such an EOS, whether we can consistently explain 
the structural properties such as mass, radius and moment of inertia as well 
as the presence of the magnetic field at the surface within the same model. 



2. Theoretical framework 



In the following, we briefly outline the procedure to obtain the nuclear 
equation of state in a non-relativistic framework and discuss its merits and 
limitations. 
2.1. Equation of state 

The phenomenological momentum and density dependent finite range in- 
teraction employed here to obtain the equation of state is a modified version 



of Seyler-Blanchard interaction [^7]]. To treat spin-polarised isospin asymmet- 
ric nuclear matter, the interaction has been generalised to include explicitly 
the spin-isospin dependent channel. The interaction between two nucleons 
with separation r and relative momentum p is given by, 



} eff 



r/a 



r/a 

(1) 



where a is the range and b defines the strength of the repulsion in the momen- 
tum dependence of the interaction. The parameters d and n are measures of 
the strength of the density dependence, and p\ and p 2 are the densities at the 
sites of the two interacting nucleons. The subscripts r and s in the strength 
parameter C TS refer to the likeness I and the unlikeness u in the isotopic spin 
and spin of the two nucleons respectively; for example, Cu refers to interac- 
tions between two neutrons or protons with parallel spins, Ci u refers to that 
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between neutrons or protons with opposite spin etc. The energy per nucleon 
E and the pressure P in the mean-field approximation can then be worked 
out[p7j as 
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Here, Jk(v) are the Fermi integrals, V® s and V 2 S are the single-particle and the 
rearrangement potentials, V^ s is the coefficient of the quadratic momentum 
dependent term in the potential and defines the effective mass m* s , T the 
temperature, and r\ is the fugacity given by r\ TS = (p TS — V® s — V 2 S )/T. For 
the unpolarised nuclear matter (NM), the expressions for V® s etc are given in 
ref . []57j . It is straightforward to extend these to the case of polarised nuclear 
matter and are as given below: 
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One usually defines the neutron and proton spin excess parameters ( spin 
asymmetry ) as 



a,, 



a.,. 



(PnT ~ Pn{)/P 

(Ppt -PpiVp • 



(5) 
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where 

P = Pn + Pp = (PnT + Pnl) + (PpT + Ppl) > ( 6 ) 

is the number density. We then define the proton fraction as x = p v j p. It is 
related to the isospin asymmetry parameter X as, 



X = (l-2x) = {p n - p p )/p . 



(7) 



We also define the spin excess parameter as Y = a n + a p and the spin- 
isospin excess parameter as Z = a n — a p . One can then express the energy 
per nucleon E/A of the NM at zero temperature as 



E/A = E V + E X X 2 + E Y Y 2 + E Z Z 2 



(8) 



where terms higher than those quadratic in X, Y and Z are neglected. Here 
Ey is the volume energy of the symmetric nuclear matter, taken as —16.1 
MeV and Ex is the usual symmetry (isospin) energy, taken to be 34.0 MeV. 
The quantities Ey and Ez are the spin and the spin-isospin symmetry ener- 
gies of the NM respectively. Their values are uncertain to some extent. We 
takepl pi E Y = 31.5 MeV and E z = 35.0 MeV in conformity with the gen- 



1/2 

eralised hydrodynamical model of UberallpO|, where {Ez/ Ex) ~ 1.1. In 
terms of the strength parameters C TS , the volume and the symmetry energies 
are written in the form, 
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(9) 



Here pp is the Fermi momentum of the one-component nuclear matter cor- 
responding to the density p of the polarised NM, A = a{(3 — 5), B = 
a{(3 - 205/9), C = a(p - 105/9), a = 8vr 2 a 3 p|/(3/i 3 ), (3 = 1- d 2 {2p) n 
and S = 6p 2 7 /{5b 2 ). For a fixed value of n, the parameters C TS , a, b and d 
are then determined by reproducing E v , E x , E Y , E z , the saturation density 
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of normal nuclear matter(po = 0.1533 fm~ 3 ), the surface energy coefficient 
(as = 18.0 MeV ), and the energy dependence of the real part of the nucleon- 
nucleus optical potential. The parameter n is determined by reproducing the 
breathing- mode energies |31"|. 



The parameters of the interaction are listed below: 

C a = -305.2 MeV a = 0.625 fm 

C lu = 902.2 MeV b = 927.5 MeV/c 

C ul = 979.4 MeV d = 0.879 fm 3n/2 

C uu = 776.2 MeV n = 1/6. 

With the above value of the parameter n, the incompressibility of symmetric 
nuclear matter is K = 240 MeV. 

2.2. Merits and limitations 

It has been tested that the above interaction reproduces quite well the 
ground state binding energies, root mean square charge radii, charge distri- 
butions and giant monopole resonance energies for a host of even-even nuclei 
ranging from 16 to very heavy systems. Interactions of this type have been 
used before with great success by Myers and Swiatecki|32|] in the context of 
nuclear mass formula. We have also seen that for symmetric nuclear mat- 
ter, our results agree extremely well with those calculated in a variational 
approach by Friedmann and Pandharipande(FP)[l2|] with + TNI inter- 
action in the density range |po < P < 2p - However, for unpolarised pure 
neutron matter, the energies calculated with our interaction are somewhat 
higher compared to the FP energies, particularly at higher densities. The 
entropy per particle for neutron matter calculated with our interaction at 
different temperatures agrees extremely well with the corresponding FP re- 
sults. In Fig.l, the energy per particle for neutron matter is displayed as a 
function of density at zero temperature. For comparision, the FP energies |fT2"| 
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and those obtained with the Bethe-Johnson(BJ) potential |33| in a sophisti- 
cated correlated basis function approach are also displayed. The BJ curve 
is very close to ours for neutron matter. This good agreement between our 



calculations and those reported in refs.fLzf and |33j suggests that the present 
interaction can be extrapolated with some confidence to neutron matter or 
to nuclear matter with large isospin asymmetry at high densities. It can also 
be mentioned that such an interaction satisfies the Landau-Migdal stability 



criteria 34 



All the well-known non-relativistic nuclear equations of state suffer from 
lack of causality at high densities. The velocity of sound in nuclear matter 
then becomes superluminal. The effective interaction used by us is no excep- 
tion. In Fig. 2, we have plotted the velocity of sound in units of c as a function 
of the ratio p/p in the case of neutron matter taking a n = 0, 0.3, and 0.5. 
It can be seen that as spin-polarisation increases, the EOS becomes softer 
and the velocity of sound v s becomes acausal only at increasingly higher den- 
sities. This superluminous behaviour of v s , particularly for the unpolarised 
neutron matter, suggests that the extrapolation of such an EOS to very high 
nuclear densities may not be advisable. 



3. Ferromagnetic phase transition 



It would be interesting to investigate whether the nuclear EOS discussed 
in the previous section, in addition to the consistent description of the nuclear 
matter and finite nuclear properties predicts a ferromagnetic transition at 
densities meaningful in the context of neutron stars. 

It has been conjectured [|35"1 , |36|1 that in the neutron star matter, in contrast 
to nuclei where the neutrons generally pair up to spin J = in their ground 
states, the neutrons may pair up to spin J = 1 at higher nuclear densities, 
thus leading to a ferromagnetic transition. To investigate this aspect, we 
calculate the energy per particle for the spin polarised neutron matter with 
a representative value of a n = 0.5 and compared with that calculated for 
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unpolarised neutron matter in Fig.l. We find that above p ~ 2.5p , the 
energy of polarised matter is lower compared to that for unpolarised neutron 
matter. This reflects that the neutrons that pair up to J = at lower 
densities undergo a transition to a spin polarised configuration as density 
builds up. In otherwords, the system prefers a ferromagnetic state for p > 
2.5 Po . 

The behaviour of the magnetic susceptibility x °f neutron matter as a 
function of density also portrays the occurence of ferromagnetic phase tran- 
sition. In general, the magnetic susceptibility is defined p7| as x — fff > where 



H is the magnetic field and M = p n {N^ — N±)/V is the total magnetisation 
per unit volume with p n being the magnetic moment of a neutron. Using the 
definition of a n [Eq.(5)], we rewrite M as M = p n a n p n . 

We need to determine the optimum value of a n using the energy minimi- 
sation criteria, 

d(E H (p,a n )/N) , 

U=a ° = °- (10) 

The total energy Eh/N per particle of a system of N number of neutrons in 
the presence of an external weak magnetic field H is, 

E H {p, a n )/N = E{p, a n )/N - {p n Ha n ) . (11) 

Expanding the energy E(p, a n ) in powers of a n upto 0(a%), we get 

x ,» T 1 2d 2 (E(p,a n )/N) . 
E(p } a n )/N = E(p,a n = 0)/N + -a 2 n 1 ^J 1 " J | Qw=0 , 

= e o + \a 2 n e 2 . (12) 

Because the energy E(p, a n ) is symmetric in a n , all the odd derivatives in 
the expansion of E(p, a n ) vanish. (It may be said here that only in the 
calculation of \i we have expanded E(p, a n ) in powers of a n , otherwise we 
have calculated it numerically.) Then, the optimum value a° is determined 
by minimising the energy [Eq. (11)] as a° = p n H/e2- Now, we can determine 
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X and it is given as, 
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Using the effective interaction given in Eq.(l), we get 
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where A\ = 1 — d 2 (2p n ) n and A 2 = PnP^/b 3 - In the limit of no interaction, e 2 
is simply given by the kinetic term alone, i.e. e^ 66 = 2 2 ^p 2 F / (3m). It is then 
straightforward to calculate the ratio Xfree/ Xi where Xfree is the magnetic 
susceptibility of the non-interacting neutron gas. The onset of a ferromag- 
netic transition is depicted by the vanishing of the ratio Xfree/ X- Further, the 
effect of the nuclear matter incompressibility K on the ferromagnetic transi- 
tion density is studied and the results are shown in Fig. 3. As the value of K 
is increased from 240 MeV to 304 MeV (by increasing the density exponent 
n of the effective interaction given by eq.(l) from 1/6 to 2/3), the density at 
which the transition takes place decreases from ~ 2.4p„ to ~ 2.3p Q . It is thus 
seen that the effective interaction given in Eq.(l) predicts a ferromagnetic 
phase transition at a density p ~ 2.4p Q . 

In previous calculations in the non-relativistic formalism using realistic 



interactions, one usually does not find such a ferromagnetic transition [18, 19 



or if there is such a possibility, it occurs at densities[19] not realisable in the 
context of neutron stars. In a relativistic framework, one finds that the sim- 



ple a + us model [20 1 does not suggest any such transition. However, with an 



improved modelpl] including other mesons like p and ir in addition to cr + cu, 
one finds a transition at about p ~ 3.5p D . But, the incompressibility of the 
normal nuclear matter obtained in this model is too high (~ 450 MeV). In 
contrast, our nuclear interaction that predicts a ferromagnetic transition at 
a similar density yields a value of incompressibility that is close to the one 
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well- accepted 0. 



4. Structural properties of neutron stars 

In this section, we explore the various static properties of neutron stars 
such as proton fraction, mass, size and moment of inertia using the proposed 
equation of state. We further study the influence of spin polarisation on 
these observables. 



4.1. Beta- equilibrium proton fraction 

In recent years, attention has been drawn to the direct URCA process 
in neutron stars which may be the primary mechanism for its rapid cooling. 
This can, however, occur only when the beta-equilibrium proton fraction x 
in the star is > 0.11, where only electrons are considered, and > 0.148, if 
both electrons and muons are considered. It would be interesting to know 
whether spin polarisation favours or disfavours direct URCA process. In 
our study, the lepton energy per particle El(p,x) is given by the relativis- 
ts, ideal Fermi-gas expression^]; in addition to e~ , p~ are also considered 



as and when they are energetically favoured. At beta-equilibrium, one has 
jj- (E(p,x) + E L (x)) = 0, where E(p,x) is the baryonic energy per particle 
including the rest masses. In Fig.4, the beta-equilibrium proton fraction thus 
obtained in the neutron star matter is displayed as a function of the baryon 
density invoking the condition of charge neutrality. The upper panel cor- 
responds to unpolarised matter, and the lower panel displays that for spin 
polarised matter with a n = 0.3 and a n = 0.5. From the upper panel, it can be 
seen that for a n = 0, x shows a peaked structure against density. When only 
e~ are considered, x increases as density increases, reaches a maximum value 
of about 0.085 at p ~ 3p Q , and then decreases to very low values at higher 
densities. With the inclusion of p~ , the structure of the curve remains almost 
the same; however, it lies higher than that of the former case, at all densities. 
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The peak value is then about 0.11. The lower panel of Fig. 4 displays pro- 
ton fractions for the spin polarised neutron star matter with a n = 0.3 and 
a n = 0.5. Here both muons and electrons are taken into consideration for 
calculating proton fractions at beta-equilibrium. With increasing spin polar- 
isation, the proton fraction becomes smaller at any density. It may be noted 
that the present EOS in use does not favour direct URCA process, since the 
proton fraction is always below the critical value. Introduction of spin polar- 
isation disfavours direct URCA process even more. It may be remarked that 
various calculations give different conclusions [1^] regarding the direct URCA 
process. It may also be mentioned that inclusion of exotic processes like 
pion condensation or kaon condensation in dense neutron star matter may 



enhance the proton fraction thereby favouring direct URCA process |39|], but 
occurence of such exotic phenomena is still very much unsettled H. 



4.2. Mass and size of neutron stars 



We now determine the structure of neutron stars using a composite EOS 
i.e. FMT, BPS, Baym-Bethe-Pethick(BBP) [|i~T| and the present interaction 
with progressively increasing densities. Then, the total mass and the size of 
the neutron star can be obtained by solving the general relativistic Tolman- 
Oppenheimer-Volkoff(TOV) equation, 



dP(r) 


G [e(r)+P(r) 


[m(r)c 2 + 4nr 3 P(r)] 


dr 




\ 2Gm(r)~ 

J- 9 

rc z 





(15) 



where, 

m(r)c 2 = [ r e(r')d s r'. (16) 
Jo 

The quantities e(r) and P{r) are the energy density and pressure at a radial 
distance r from the centre, and are given by the equation of state. The mass 
of the star contained within a distance r is given by m(r). The size of the star 
is determined by the boundary condition P(R) = and the total mass M 
of the star integrated upto the surface R is given by M = m(R). The single 
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integration constant needed to solve the TOV equation is P c , the pressure at 
the center of the star calculated at a given central density p c . 

The mass functions of the star thus obtained as a function of its cen- 
tral density are shown in Fig. 5 for four different values of spin polarisation 
(a n = 0.0, 0.3, 0.4, 0.5). The radii, central densities and surface redshifts 
corresponding to the maximum mass M max configuration are tabulated in 



Table 1 for three values of a n . The surface redshift z s is defined JTT| as 

2GMr^_ L 



Rc 2 

Values of M max /M & are also given in the same table. It can be seen that 
the maximum mass and the corresponding radius and surface redshift z s 
decrease with increasing polarisation. On the other hand, the central density 
pertaining to M max configuration increases as a n increases. This is because of 
the fact that the spin polarised neutron star matter is more compressible than 
the unpolarised one. The measured mass [[40H of 4£/0900 — 40, (1.85±O.3)M 
possibly provides the lower limit of the maximum mass of the neutron star. 
Our calculations with the values of a n taken are well within this limit. The 
maximum mass of neutron star obtained from calculations with a n = 0.6 
and above are found to be below the present observational limits and hence 
we have restricted our calculations upto a n = 0.5. It may be mentioned that 
the unpolarised neutron matter is close to being superluminal at the central 
density corresponding to the maximum mass. For polarised neutron star 
matter though the central density increases significantly, the sound velocity 
is always luminal because of the softness of the polarised matter towards 
compression. 

We have also studied the sensitivity of the mass distribution m(r) to a n , 
where m(r) denotes the total mass contained within a given radial distance r 
from the centre. Fig. 6 displays m(r)/M max as a function of the density p(r) 
at that point r, where M max corresponds to the maximum mass configuration 
at a given ot n . The zero of the abscissa refers to the surface of the neutron 
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star; the points where the mass functions m(r) meet the abscissa refer to the 
centre of the star. It is found that nearly 90% of the mass of the neutron 
star lies at densities higher than 3p a . 



4.3. Moment of inertia of neutron stars 

The moment of inertia of neutron stars is calculated by assuming the star 
to be rotating slowly with an uniform angular velocity f2 ||41|| . The angular 



velocity u(r) of a point in the star measured with respect to the angular 
velocity of the local inertial frame is determined by the equation, 



]_d_ 

r 4 dr 



4 .duj 

r J-r- 

dr 



+ -fo = 0, (18) 
r dr 



where, 

1/2 



j = e~^ l- 2Gm 9 {r) ) . (19) 



rc 2 



The function 0(r) is constrained by the condition, 



/ 2GM 

(r V(r) = constant = fi(R) \j 1 - (20) 



where the chemical potential fi(r) is defined as, 



= ! m + p( i ) 

p(r) 

Using these relations, Eq.(18) can be solved subject to the boundary condi- 
tions that u{r) is regular as r — > 0, and u{r) — > Q as r — > oo. Then 
moment of inertia of the star can be calculated using the definition / = J/Q, 
where the total angular momentum J is given as 

Values of / thus obtained for three values of a n are plotted as a function 
of the central density in Fig. 7. It can be seen that the maximum value of / 
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greatly depends on a n . As a n is increased from 0.0 to 0.5, I max has decreased 
by about 50%. 

We have thus calculated the important structural properties of neutron 
stars and studied their dependence on the spin polarisation factor. It is 
found that upto a value a n = 0.5, the properties obtained here are well 
within the observational limits. Motivated by this success in explaining 
the structural properties of neutron stars, we are curious to see how far 
we can account for the surface magnetic field. Our present study suggests 
that the neutron star matter at densities p > 2.5p Q is spin polarised. Taking 
a n = 0.5, one immediately realises that the surface magnetic field is largely 
overestimated(~ 10 16 G). As noted earlier, it is probable that that at core 
densities one finds quark matter rather than spin polarised nuclear matter. 
Therefore, we in the following section explore this plausibility and its impli- 
cations upon both the structural and magnetic properties of stars. 



5. Hybrid stars 

Here, we construct the (3— equilibrated, electrically neutral quark matter 
equation of state and employ it to understand the structural properties of 
the hybrid stars,i.e. quark cores with nucleon envelopes. 

The equation of state of a three flavour quark matter consisting u, d and 
s quarks is obtained here using the phenomenological MIT bag model [f!2]|. 
The total kinetic energy density of a system of non-interacting, relativistic 
quarks of flavour r and mass m T is given as, 



3 (m T c 2 ) 



7T (He) 



x T Jl + x 2 T (1 + 2<) - ln(x T + Jl + x 2 T 



(23) 



where x T = p T F /(m T c), p T F being the Fermi momentum and is related to the 

1 /3 

quark number density p T of a given flavour as p T F = h(n 2 p T ) . The densities 
pertaining to the three flavours can be expressed in terms of the total quark 
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number density p q and the asymmetry parameters S u d and S us as: 

Pu = (Pg/3) [1 - 5 ud - 5 US ] , 
Pd = (Pg/3) [1 + 2S ud - S us \ , 

p s = ( Pq /3)[l-5 ud + 25 us ], (24) 

where 5 ud = (pa - Pu)/p g , S us = (p s - p u )/p q and p q = p u + p d + p s . Similarly, 
the energy density pertaining to a system of relativistic non- interacting 
electron gas can be calculated^ . 

We then determine the equilibrium composition of the quark matter sub- 
ject to the (3— equilibrium condition, 

Pd~ Pu = Pe and p d = p s , (25) 
and the charge neutrality condition, 

Pe = - (2p u - Pd- Ps) ■ (26) 

Using Eq.(24) one obtains, p e = — (p g /3) (S u d+S us ). Similarly, we can express 
the chemical potentials p u , pd and p s in terms of the three quantities p q , 5 u d 
and (Lo as follows: 



9e q \ 


de q 


1 + 5 ud de q 


1 + 8 US de q 


dpu) 


\p d .p. Q pq 


Pq d5 ud 


P q dS us 


deA 


= d_e ± _ 


1 - 5 ud de q 


S us de q 


dpd) 


\ P u,Ps 


P q d5 ud 


P q d5 us ' 




_ de q 


Kd de q ! 1 


- S us de q 


dps) 


lpu ' Pd ~ dp q 


p q d5 ud 


p q d5 us ' 



Pu 

Pd 



^^q _ ^^q _ ^ua ^'-q "us ^>-q ^J) 

^ \d P s)^ Pd d Pq Pq dSud p q d5 u ; 1 ; 

where e q = J2 T e r + B is the total quark energy density and B is the bag 
parameter. Using these expressions, the (3— equilibrium conditions can be 
rewritten as, 

de i de q = Q 



dKd d5 us 

deq 4. 1 deq - u (?R) 



Pq dSud Pq dS v 
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where p e = ypp^c 2 + rn^c 4 . Thus, for a given baryon density p& = p g /3, 
the three quantities p e , 5 u d and 6 US are fixed by the Eqs.(26) and (28). Sub- 
sequently, the equation of state is completely described by the total energy 
density €qm and the pressure Pqm of the system calculated for a given p b 
using the definitions, 

€ QM = ^ e r{Pq, Sudi 8 U s) + e h{Pqi ^udi $us) + B, 

T 

Pqm = p q e QM . (29) 

It may be mentioned that in our present study we have omitted the lowest 
order quark-quark interaction terms as it can be effectively absorbed into the 
bag constant f43|. The masses of the quarks are taken to be: m u = 5 MeV, 
m d = 10 MeV and m s = 200 MeV. 

To know whether there is a phase transition from the nuclear matter to 
quark matter, we compare the total energies per baryon obtained in the two 
phases. The EOS corresponding to NM is given by Eqs.(2-3) and that of 
QM is given by Eq.(29). In Figs. 8 and 9 we show the energy per baryon as 
a function of the baryon density for two values of bag constant, B 1 / 4 = 155 
MeV and B 1 / 4 = 170 MeV respectively. These are compared with the curves 
obtained using a n = 0.0 and a n = 0.5 in the same figures. The baryon 
density p& at which the polarised(a„ = 0.5) and an unpolarised(a n = 0.0) 
curves intersect is denoted by Pfm- For the nuclear interaction given in 
Eq.(l), pfm — 0.405 fm~ 3 . Similarly, the baryon density at which either of 
the NM curves intersect with the QM one is denoted by Phq- 

It can be seen in Fig. 8 that for B 1 ^ = 155 MeV, the unpolarised NM is 
energetically favoured upto a density p& ~ 0.255 fm~ 3 . As p b is increased 
further, the quark matter is found to be the lowest energy state. Therefore, 
there is no region of spin polarised NM in this particular case. On the other 
hand, for 5 1 / 4 = 170 MeV, it can be seen that for densities within the values 
p FM — 0.405 fm~ 3 and puq — 1-05 fm~ 3 , the polarised NM is the state 
of lowest energy. Thus, in this case(Fig.9), one has a spin polarised region 
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sandwitched between a quark matter core (p^ > Phq) and a unpolarised nu- 
clear matter envelope^ < Pfm)- Further, we also studied the dependence 
of the width of the spin polarised region, roughly given by (phq — Pfm), 
on the strange quark mass. In Fig.10, we have plotted (phq — Pfm) I Phq 
as a function of the bag paramter for three values of m s . Negative values of 
{.Phq— Pfm) I Phq indicate that there is no region of spin polarised NM inside 
a hybrid star. In otherwords, to have a spin polarised region sandwitched 
between an unpolarised NM envelope and a QM core, p HQ must be greater 
than pfm- It can be seen that for allowed values of B l / A and m s , the presence 
and non-presence of a polarised region are equally probable. Thus, it is clear 
from the above discussions that the allowed range of B and m s values is not 
able to decide upon the presence/non-presence of a spin polarised region in 
the hybrid star. It would be interesting to know whether the observational 
limits on the mass and the size of pulsars would help in finding an answer. 

In view of this, we explore the structural properties of the hybrid stars 
choosing three particular values of the bag parameter B while keeping the 
strange quark mass m s fixed at 200 MeV. The three different configurations 
corresponding to the three sets of MIT bag model parameters are shown 
schematically in Fig.ll. It is clear from the figure that for a fixed value 
of m s , one can appropriately choose the value of B so that the calculated 
surface magnetic field is of the order of the observed value(~ 10 12 G). The 
mass, size, central density and surface redshift were then calculated for the 
three values of B. The results obtained pertaining to the maximum mass 
configuration are given in Table 2. It can be seen that M max /M Q and R 
decreases as the value of B is increased. For B 1 ^ = 170 MeV, M max /M Q is 
found to be less than the observed value. Therefore, one may conclude that 
large regions of spin polarised matter in a hybrid star is ruled out by the 
observational limit on the mass of pulsars. But, our present study cannot 
decide upon either of the two configurations obtained with B 1 ^ = 155 MeV 
and I? 1 / 4 = 164 MeV, since the mass of the respective stars and their radii 
are well within the acceptable limits. We would however like to mention that 
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using the proposed equation of state and the MIT bag model with acceptable 
parameters^ 1 / 4 ~ 164 MeV, m s =200 MeV), one can in principle describe 
both the structural properties as well as the surface magnetic field satisfac- 
torily. 

6. Summary 

To summarise, we have constructed a nuclear equation of state from a 
finite range momentum and density dependent interaction and then applied it 
to investigate some properties of spin polarised nuclear matter with particular 
reference to the neutron star matter. The parameters of the interaction have 
a firm basis in the well known properties of nuclear matter and of finite nuclei. 
Extrapolating this interaction to neutron matter and to higher densities, it 
is found that the present equation of state agrees well with those obtained 
from more sophisticated calculations. 

Introducing spin degrees of freedom, it is seen that at density p ~ 2.5p Q , 
the neutron star matter undergoes a ferromagnetic phase transition. This 
aspect is demonstrated by studying the density behaviour of the total en- 
ergy(Fig.l) and the magnetic susceptibility (Fig. 3). To the best of our knowl- 
edge, this is the first non-relativistic calculation that gives a ferromagnetic 
transition at a density well realisable in the neutron star core. 

The proposed equation of state is then applied to the investigation of the 
structure of neutron stars. With increasing polarisation, the maximum mass 
and the corresponding radius decreases whereas the value of central density 
increases. The maximum mass of stars and their radii obtained from calcula- 
tions with a n < 0.5 are found to be well in agreement with the observations. 

There is a good possibility that one finds quark matter rather than spin 
polarised nuclear matter at the core of stars. We therefore investigated this 
plausibility using the MIT bag model. We found that using the proposed 
EOS and the MIT bag model, one can in principle obtain the maximum 
mass, size and the surface magnetic field of hybrid stars well within the 
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acceptable limits. 
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FIGURE CAPTIONS 



Fig.l. The energy per particle for pure neutron matter as a function of den- 
sity is plotted for a n = and a n = 0.5. The calculated results of Friedman- 



Pandharipande[12] and with Bethe- Johnson potential [33] are also shown 



Fig. 2. Velocity of sound v s obtained in units of c for pure neutron matter 
taking a n = 0.0, 0.3 and 0.5 is plotted as a function of the density ratio pj p a , 
where p Q = 0.1533 fm -3 is the symmetric nuclear matter saturation density. 

Fig. 3. Magnetic susceptibility x °f pure neutron matter calculated for two 
values of nuclear matter incompressibility K is shown as a function of the 
density p, where Xfree is the magnetic susceptibility of non-interacting neu- 
tron gas. 

Fig. 4. Beta-equilibrium proton fraction x as a function of density of neutron 
star matter. In the upper panel, results are shown for a n = with e~ and 
e~ + pT considered for f3— equilibrium. In the lower panel, the results for 
a n = 0.3 and a n = 0.5 are shown with + pT . 

Fig. 5. The neutron star mass is plotted as a function of central density. 

Fig. 6. The integrated mass upto radius r is plotted as a function of the 
density at radius r for different spin polarisation a n . 

Fig. 7 Moment of inertia obtained using a ra =0.0, 0.3 and 0.5 is shown as 
a function of density p, where p Q = 0.1533 fm -3 is the symmetric nuclear 
matter saturation density. 

Fig. 8 The total energy per baryon of quark matter(QM) calculated us- 
ing the bag model picture is compared with the energies per baryon of 
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unpolarised(a n = 0.0) and polarised(o; n = 0.5) nuclear matter(NM). 
Fig.9 Same as Fig.8, but for S 1 / 4 =170 MeV. 

Fig. 10 The width of the spin polarised region, roughly given by (phq — 
Pfm)I Phq, is plotted as a function of the bag parameter B for three values 
of strange quark mass, m s =100, 200 and 300 MeV. Negative values of (phq — 
Pfm)I Phq imply that there is no spin polarised region inside a star. 

Fig. 11 Schematic representation of the three configurations considered in 
our study of hybrid stars. 

TABLE CAPTIONS 

Table 1 Values of the mass M max , size R, central density p c , surface redshift 
z s and moment of inertia / obtained in the case of neutron stars using three 
values of a n are shown. The tabulated values correspond to the maximum 
mass configuration. 

Table 2 Values of the mass M maxi size R, central density p c , surface redshift 
z s and moment of inertia / obtained in the case of hybrid stars using three 
values of bag parameter B are shown. The strange quark mass m s =200 MeV. 
The tabulated values correspond to the maximum mass configuration. 
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Table 1 



OLn 


M max /M e 


R 

(km) 


Pc/Po 


Z s 


0.0 


2.03 


10.3 


7.6 


0.55 


0.3 


1.89 


9.9 


8.4 


0.52 


0.5 


1.53 


8.7 


11.7 


0.44 



Table 2 



£1/4 


M max /M Q 


R 


Pc/Po 




(MeV) 




(km) 






155 


1.62 


10.1 


8.8 


0.38 


164 


1.46 


9.3 


10.3 


0.37 


170 


1.37 


8.8 


11.4 


0.36 
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